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1 Introduction 

The last ten years have witnessed tremendous developments in the control of 
systems at the nanoscale, allowing the fabrication of artificial structures that are 
effectively two-dimensional [l|; one-dimensional [1; 0] or zero-dimensional 
Simultaneously, progress in optical trapping of atoms have also lead to the pos- 
sibility of realizing quantum systems of controllable effective dimensionality Q. 
From the theoretical point of view, these low dimensional systems are expected to 
exhibit a rich physics resulting from quantum coherence which is now becoming 
experimentally accessible. A particularly striking property in low dimension is the 
rS ] phenomenon of Anderson localisation 0] , in which disorder makes the wavefunction 
of particles localised. In both one Q and two Q dimensions, any small amount of 
disorder is expected to localise all particles. In one dimension, random matrix the- 
ories have provided a highly detailed understanding of the localisation phenomenon 
[Iol | , while in two dimensions the regime of weak localisation has been thoroughly 
studied. Our group has been exploring at a theoretical level the physics of local- 
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isation in low-dimensional artificial structures. Two main lines of research have 
been emphasised. In the first one, HI; HI, reviewed in Sec.[2j mesoscopic samples 
of metallic spin glasses are analysed. In these systems, the underlying disorder is 
generated by the spin glass state, and weak localisation is used as a tool to probe 
its statistical properties. The fluctuations of observable quantities such as the con- 
ductance and the correlations between them in the weak localisation regime are 
shown to provide information on the overlap between the states of the spin glass 
of the spin glass state. In the second line of research, reviewed in Sec. [3j trapped 



bosonic atomic gases are considered. |11|; |12|; [15| In these systems, in contrast to 



the spin glass systems, there is a high degree of control of the disorder acting on 
the system, and the aim is to explore how Anderson localisation is affected by the 
interaction and the statistics of the quantum particles, a long standing problem of 
condensed matter physics. [l6| 



2 Electronic transport through a spin glass nanowire 

2. 1 Short overview of spin glass physics 

Spin glasses are amorphous magnetic phases. The canonical examples of such 
glasses are obtained by doping with magnetic impurities noble non-magnetic met- 
als, such as Cu:Mn, Ag:Mn, Au:Fe. At high enough concentrations, the coupling 
between the randomly located impurities lead to the appearance of a new mag- 
netic phase. In this phase, the spins are frozen but lack any simple long range 
order. The analogy with conventional glasses is at the origin of the denomination 
of this magnetic phase. The understanding of the nature of this frozen phase and 
the associated transition from the paramagnetic phase has motivated numerous 
experimental, theoretical and numerical works over the last three decades. Ex- 
perimentally, this freezing manifest itself through the sharp reduction of magnetic 
susceptibility x'( w ~ * 0,T) below the glass transition (but not on e.g the specific 
heat). This decrease of x'( w — > 0) i s accompanied by a broad spectrum of the 
imaginary part of manifesting the very slow magnetic relaxation in the spin 

glass phase. Indeed, no intrinsic time scale characteristic of this relaxation can be 
determined experimentally : the spin glass is said to be aging. Its only relevant 
time scale is its age t w : correlation and response functions are not stationary and 
depend explicitly on both times t w and t w + t. This dependence on history of the 
spin glass properties possesses remarkable properties : for example, the spin suscep- 
tibility ~ * 0,T) crucially depends on the presence or not of a small magnetic 
field upon cooling through the transition temperature (field cooled or zero field 
cooled protocols). The response of the aging properties of the thermo-remanent 
magnetisation to small variations of temperatures manifest the so-called memory 
and rejuvenation effects. 

These unusual properties have focused the attention of a large community of 
condensed matter theoreticians. Most studies have focused on the Ising spin glass, 
expected to describe the low energy excitations in the presence of spin anisotropy 
(originating from e.g. spin-orbit coupling). For Ising spin glasses, G. Parisi pro- 
posed a mean field Ansatz which is now understood to be exact 17; 18; 0|. In 



this mean field solution for the statics, the phase space is complex with different 
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crgodic components (see [2CJ; [21| for pedagogical reviews). These crgodic com- 
ponents correspond to different pure states. Naively, we expect from such a mean 
field scenario that whenever we enter the spin glass phase by crossing the transition 
temperature, we end up in a randomly chosen ergodic component where relaxation 
dynamics takes place. If we repeat the quenching procedure in the same sample, 
we can end up in another of these ergodic component. History dependence fol- 
lows directly. The number of different ergodic components (or "energy valleys") 
depends of the complexity of the phase space, and is encoded in the distribution 
of probability of distances between these pure states. These distances are conve- 
niently encoded using the overlap between spin configurations. Let us consider two 
classical spin configurations }j and {s!> }i (same positions of the spins, but 
different orientations). The overlap between these configurations is defined as 

. N tmp 

(1) <2l2 = jt— Si^-S^, 



where Ni mp is the number of magnetic impurities (spins). For Ising spins, this 

overlap will count by how many spins flips the configurations {S^- }j and {S^ }j 
differ. The distribution of overlaps is a central object in the Parisi's mean field 



solution [20|; [21| , and its distribution is indeed the proposed order parameter for the 
spin glass transition. Moreover besides allowing to encode part of the complexity 
of the phase space, this overlap is also a useful tool to measure variations of a spin 
configuration under variation of an external parameter (temperature or magnetic 
field), or as time evolves, for which configurations 1 and 2 corresponds to different 
times (see e.g (33|). 

Coming back to the mean field solution, its relevance to the physical situation 
D = 3 is still an active domain of research. Besides this static solution, solutions 
of mean-field dynamics have been obtained by Cugliandolo and Kurchan (see [33| 
for a recent detailed review). Mean field aging dynamics can also be described 
in the so-called wcak-crgodicity scenario proposed by Bouchaud [28[. Moreover, 
even the complexity of the phase space have been questioned. In another proposal 
based on scaling ideas 0; [4(| , the peculiar properties of spin glasses are attributed 
solely to the aging dynamics, without any need for a complex phase space. In 
this phenomcnological theory, the slow dynamics is described analogously to the 
evolution of a quenched ferromagnet : domains of various sizes and shapes evolve 
slowly in time. It is fair to say that, while these numerous results have deepened our 
understanding of the amazing properties of spin glasses, a complete understanding 
of real spin glasses is still lacking. The purpose of our work is to propose another 
experimentally accessible probe of metallic spin glasses, allowing interesting and 
crucial comparisons with current theories. 



2. 2 What does coherent electronic transport probe in a spin glass ? 

In this work, we take advantage of the metallic nature of canonical spin glasses 
and consider their electronic transport properties in the coherent regime. In metals, 
electron's inelastic scattering off e.g. phonons limit the phase coherence of these 
electrons. This dephasing can be described phcnomcnologically by a dephasing 
length scale L^T) beyond which phase coherence can be safely neglected when 
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describing transport properties. At room temperatures, the high density of phonons 
induce an L<f,{T) much smaller than conductors sizes : the classical Drude theory 
correctly describes the conductance of metallic wires. However upon lowering the 
temperature, L^iT) increases up to a few /im, allowing to study transport in wires 
of size comparable to L$(T). When describing the electronic transport through 
such devices, interferences effects have to be taken into account : this is the so- 
called mesoscopic regime (for a recent and very detailed introduction to this field, 
see 113]). In this regime, the conductance is sample dependent. Moreover, when 
applying a small magnetic field transverse to the wire, the conductance fluctuates 
in a random manner. The corresponding magneto-conductance trace is considered 
as a fingerprint of the configuration of disorder in the sample : two samples with the 
same amount of disorder differ by the exact realization of the random potential. 
Hence they display magneto-conductance traces that are different, but with the 
same amplitude. 

In a metallic spin glass, part of the randomness consists in the random orien- 
tations of the impurities magnetic moments. The influence of these random spins 
on the fluctuations of the conductance was soon realised by Al'tshuler and Spivak 
[25| and specific consequences of the droplet theory of spin glasses were discussed 
in [371 ]. Experimentally, dc Vcgvar et al. pioneered the mesoscopic transport stud- 
ies in metallic spin glasses by clearly demonstrating the existence of reproducible 
magneto-conductance traces below the glass transition in Cu:Mn wires [3^|. By 
using the antisymmetric component of conductance upon time-reversal symmetry, 
they attributed the appearance of these coherent conductance fluctuations to the 
freezing of impurities spins. In parallel, Israeloff et al. focused on the electrical noise 
in s pin glas s Cu:Mn wires, i.e. on the time dependence of conductance fluctuations 
[24l; |44|; 15^1 . The dynamical characteristics of these spontaneous fluctuations were 
analysed in view of both the mean-field and the droplets theories of spin glasses 
(see [H^ for a review). However, comparison with existing theoretical ideas 2^ 37 1 
was hampered by the necessity of averaging over the conductance distribution in 
theoretical predictions. Later, universal conductance fluctuations were analysed in 
diluted magnetic semiconductor in the spin glass regime |47j |. In these compound, a 
rapid increase of the amplitude of these fluctuations was noted below the spin glass 
transition temperature, as well as the usual field cooled/zero field cooled depen- 
dence of these fluctuations. More recently, similar studies have been conducted on 
Au:Fe wires [IH Hij. These different studies have validated the use of mesoscopic 
conductance fluctuations as an interesting new tool to probe the spin glass freez- 
ing. Most of this work have focused on either the amplitude of these fluctuations, 
or their time-dependence (electrical noise). In the present study, we consider the 
correlation between traces of magneto-conductance as a new and most relevant way 
to analyse spin glass properties. 



2. 3 Probing the distribution of overlaps through conductance correlations 

We now focus on the description of electron's coherent transport in a metallic 
spin glass. In this paper we will focus on the spin glass phase, and motivated 
by previous experimental results we will assume that the dephasing length scale 
L<f,(T) (which is dominated by spins able to flip during the diffusion time scale) is 
comparable with sample size. The spins that do not contribute to L^(T) are then 
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described by frozen classical spins with random orientations (we do not assume 
these spins to be Ising like). We consider lattice model describing the electrons in 
the presence of these impurities which contribute both to a scalar potential m and 
a magnetic part Si : 

(2) H = ]T t 4 tS c jlS + ^Vi c\ s a, s +Y1 J ?ss ' -^ C L C M' 

The transport properties of this model in the coherent regime were analysed both us- 
ing diagrammatic perturbation theory 14 ; 2|| , and a numerical Landauer approach 



1 1 31 ] - This last technique consists in obtaining the conductance of a finite sample 
for a given configuration of disorder, i.e without any averaging procedure. Follow- 
ing Landauer formula, this conductance is deduced from the electron's scattering 
matrix of the system connected to two reservoirs. In this numerical procedure, the 
scalar site-disorder Vi are uniformly distributed in Vi £ [-W/2, +W/2] in units of 
t = 1. We typically used W of order 0.5. The spin glass configuration is mimicked 
numerically by choosing spins Si randomly on the unit sphere, without any spatial 
correlations. As explained below special care was devoted to correlations (overlap) 
between different spin configurations. 

Let us first review the situation without magnetic disorder (J = 0). For a 
nanowire of size L yi L z -C ~ L Xl coherent electronic diffusion is effectively 
one-dimensional (ID). In this quasi-lD regime, the Anderson's localisation length 
scales proportionally to the number of propagating modes N e : £ = N e l e where l e 
is the elastic mean free path. For wire's length L x much smaller than £, in the so- 
called weak-localisation regime, the distribution of conductance P(G) is Gaussian. 
Without magnetic impurities (J = 0), its variance reads ((8G) 2 ) V = s 2 F(L/L l p) 
with [H (see Fig. HJ 

(3) F(x) = 12(-\ E m + X 2 } 2 - -i (* 2 csch2x + xcoth(x) 2) , 

where ()y denotes an average over the random scattering potential V(r) = ViS(r— 
Yi) , and SG is defined relative to the average conductance SG = G — (G)v- In this 
formula s is the level degeneracy, which corresponds to s = 2 for spins 1/2 . For 
an infinite dephasing length scale (x = 0), we recover the universal conductance 
value ((5G) 2 ) V = s 2 ^, whereas ((<5G) 2 ) y vanishes as expected in the classical 
limit L ^> L^. 

In the presence of magnetic impurities, the conductance G is now a function of 
both the configuration of scalar potential V(r) = J^i v i^( r ~ r i) an d the configura- 
tion of frozen spins {Si}. For a given configuration of spin {Si}, we can consider 
the distribution of G(V, {Si}) as V varies. This distribution Py(G({Si})) is still 
gaussian and independent of the configuration of spins. The length dependence of 
its variance is modified by the presence of the symmetry breaking disorder. A new 
magnetic dephasing length L rn appears (here given to order J" 1 ) : 

(4) L 2 m = hD(Ann lmpP oJ 2 S 2 y l , 

where D is the diffusion coefficient in the sample, nj mp the density of impurities 
and po the electron's density of states. In the limit L mi L L<p, the variance of 
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Figure 1 Function F(x = L/L$) describing the scaling of the variance of the con- 
ductance distribution without magnetic impurities : see eq. ((3]). The function H(x) 
represented on the right of the figure, describes the scaling of this variance when magnetic 
impurities are present, and a new dephasing length L m appears : see eq. ©. 




x=L/L<|) x=L/Lm 



P V {G) reads 

([5G{{S l })f) v = H{L/L m ) (5) 

= If (o) + -f (*±) + i F (V2—) + -f (J*±) . 

4 ivSW 4 V W 4 \\3F m ) 

This function H(x) is represented in figure [T] for — > oo. 

We now focus on the dependence of this conductance on variations of the con- 
figurations of spins. For a given configuration of V, we consider the variation of 
conductance between two spin configurations (which differ by the orientation of the 
frozen spins {Si}, but not their position) : 

(6) G(y,{s^})-G{v,{sf }). 

As V is varied, this quantity fluctuates similarly to G. The corresponding distri- 
bution, obtained using the numerical Landauer method, is represented in Fig. [2] 
for several pairs of configurations {Sj } and {Sj }. As expected, this distribution 
is still gaussian, with a vanishing average. For a fixed geometry, its variance ap- 
pears to be parameterised by the overlap (pT|) between the two spin configurations : 
the distributions of G for different pairs of configurations with similar overlap col- 
lapse. Indeed, its variance is entirely encoded by the average correlation between 
the conductances in two different frozen spin configurations , {S^ } and {S^} : 

(7) (AG)| (1)S(2) = (6G (V, {S^}) SG (y, {^f } }) ) y , 

where SG(V, {Sf } }) = G(V, {Sf '}) - (G(V, {Stilly- The length dependence of 
these correlations is obtained by extending ([5]) [lj; [25[ . In the same limit L m <C Lj, 
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Figure 2 Distribution of the variation of conductance between two spin configurations 
{Sj } and {Sj } in a given sample V. These distributions are obtained using a numerical 
Landauer approach for the model ([2]) . The coupling J between the spins of the electron 
and the impurities was chosen as J/t = 0.1, the system size was 40 x 1600. For each config- 
uration of spins, 5000 configurations of scalar disorder V were generated. The histograms 
over these 5000 points are shown for 5 pairs of spin configurations. While the distributions 
for pairs with a spin overlap Q ~ 0.91 collapse within a good approximation, the distri- 
bution corresponding to Q = 0.99 clearly deviates from these other distributions. This 

result illustrates the parameterisation of the distribution of G fv, {S^jj — G (v, {S^}J 
by the spin overlap Q. 

J=0.1,L = 1600, Nv = 5000 points 



♦ Pair 1, Q - 0.9899 
Gaussian PDF 

♦ Pair 2, Q = 0.9127 

♦ Pair 3, Q = 0.9100 

♦ Pair 4, Q = 0.9078 

♦ Pair 5, Q = 0.9098 
Gaussian PDF 
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as in ([5]) , it reads 

(A«)l,. s <» = ^(^e£) + ^(^|£) 

where Q is the overlap between the two spin configurations defined in ([T]). This 
result is the starting point for our proposal to measure the overlap between random 
spin configurations : in a given sample when orientations of the spins are varied, 
the variations of the quantity © arc entirely encoded by the overlap between the 
corresponding spin configuration. Such a measurement would allow for the first 
experimental determination of spin configuration's overlap. 

A crucial step needed to experimentally measure ((6]) is to sample the corre- 
sponding distribution. This distribution was defined by considering the statistics 
of G when the potential V is varied. However, one cannot vary the diffusing po- 
tential in a semiconducting spin glass without modifying completely the spin glass 
configuration ! We thus have to resort to the usual ergodic hypothesis of weak lo- 
calisation (see [fiH for a recent discussion). Following this hypothesis, the statistics 
of G are usually determined experimentally by considering a magneto-conductance 
trace of a given metallic sample. When a flux is applied perpendicular to the sam- 



8 D. Carpentier, E. Orignac, G. Paulin and T. Roscilde 



pie, the conductance fluctuates in a random but reproducible way (see Fig. - The 
variance of the corresponding distribution approximates correctly the variance of 
the distribution P V (G) [H. 

Figure 3 Magnetoconductance traces obtained by a numerical Landauer approach 
(la ). The different traces correspond to the same sample, i.e configuration of scalar disor- 
der, but different configurations (i.e orientations) of the frozen spins. 




Flux Quanta Through the Sample 

Figure 4 Proposed protocols to measure the correlations (overlap) between spin con- 
figurations corresponding to different quenches or different waiting times. In the first case, 
the different spins states {S[ ■ } are labelled by the index of the quench n. All these states 
correspond to the spin configuration obtained after a time t w after entering the spin glass 
phase. At this time t w , the spin configuration is sampled through the measurement of 
the corresponding magnetoconductance trace. The correlation between these traces gives 
access to the overlap between the corresponding spin states. In the second case, the spin 



configuration correspond to different waiting time t 
through magnetoconductance measurements. 



Similarly, these states are probed 



Temperature, T 



different spin config. 
= different quenches 

{sf +2) l 



Temperature, T 




different spin config. 
= different waiting times 

-(n+2). 




n + 2 



n + 2 



Following this line of reasoning, the proposed quantity to measure instead of ||7J) 



is 



(SG(y,{S^})SG(y,{sf})) B 

1 



5G(V, {S B)SG(V, {S) 2) }, B)dB (9) 
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A crucial approximation made in using as an evaluation of ([7]) is the neglect of 
the spin glass response to the (weak) magnetic field. A crucial step in this possible 
experimental route would thus consist in determining the optimal maximum field 
-B max below which this spin glass response can be neglected, but high enough to 
provide the necessary sampling of the statistics of G. 

In conclusion, we have proposed a first route to measurements of spin configura- 
tions overlap in metallic mesoscopic conductors. This route consist in determining 
the magneto-conductance traces for different spin configurations. The correlation 
between these traces gives access through eq. © to the overlap between the cor- 
responding spin configurations. In a spin glass, these different spin configurations 
can be accessed either by successive cooling below the glass transition, or e.g by dif- 
ferent waiting times. The proposed protocols corresponding to these two situations 
are illustrated in Fig. [4j 



3 Anderson localisation with cold atoms 



3. 1 Introduction 

The field of ultracold atoms has emerged in recent years as a new research 
field at the frontier between atomic, molecular and optical physics and condensed 
matter physics. @; 0; |49| . Trapping potentials of arbitrary shapes can be Fourier- 
synthesised via the superposition of different laser standing waves. In particular 
a single standing wave confines the atoms into planes, with negligible inter-plane 
coupling (over the time duration of the experiments, typically lOOms-ls); two or- 
thogonal standing waves confine the atoms into tubes, while three standing waves 
create a cubic optical lattice. More sophisticated lattice geometries can be then 
achieved by superposing further spatial Fourier components. Hence optical poten- 
tials give unprecedented access to low-dimensional quantum systems and to arbi- 
trary lattice Hamiltonians. When ultracold atoms are loaded in the lowest band of 
optical lattices, their physics is faithfully described by Bose- and Fermi-Hubbard 
models. The depth of the optical lattice allows to control the inter-site hopping, 
and hence the relevance of interactions against the kinetic energy over several or- 
ders of magnitude. Interactions are generally short-ranged (on-site only), and can 
be further tuned with the use of Fcshbach resonances [59( . Hence experiments can 
span all regimes of lattice Hamiltonians, from the weakly interacting one to the 
strongly interacting one, and probe fundamental phenomena such as the onset of 



Mott insulating physics and the suppression of quantum coherence [42l ; [61 1 . 

Yet cold atomic samples contain typically 10 5 — 10 6 atoms; when standing waves 
are applied to divide the system into planes or tubes, each layer/tube contains typ- 
ically 10 3 — I0 4 / 10-100 atoms respectively Hence trapped cold atoms can be 
viewed as a particular instance of mesoscopic systems (43| , where the finite system 
size, dictated by a global (parabolic) trapping imposed on the atoms, can be far 
exceeded by the quantum coherence length. Coherent quantum many-body phe- 
nomena out of equilibrium can be probed in the system via a real-time control 
on the Hamiltonian parameters, and they have been demonstrated in recent ex- 
periments on Bloch oscillations in periodic lattices [3], Josephson oscillations in 
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a double- well potential [23[, coherent dipole oscillations in a parabolic trap [381 ]. 
transitions from integrability to non-integrability etc. 

Our recent theoretical activities on mcsoscopic aspects of trapped cold atoms 
have been focused on the physics of Anderson localisation Q of interacting bosonic 
atoms. Anderson localisation of interacting bosons is a long-standing problem in 
condensed matter physics [lit > an d various experimental realizations have been at- 
tempted with 4 He in a porous medium such as Vycor, aerogels or xerogels [3(| or 
with disordered Josephson junctions j65|. However an experimental test of funda- 
mental theoretical predictions is still lacking. Thanks to the flexibility in designing 
the trapping potential, cold atomic systems offer the possibility to explore one di- 
mensional geometries, in which strong localisation effects are dramatic [1| and have 
been extensively studied theoretically [H(J. Moreover, the use of optical lattices 
and Fcshbach resonances allow to study the interplay of disorder and a continu- 
ously tuned interaction. 

Two methods to realise a controlled random (or pseudo-random) optical poten- 
tial are currently available. One consists of imaging a disordered diffusive plate onto 
the atomic sample, thus creating a so-called laser-speckle potential [3ll; [32; [ill; IHlj . 
The reduction of the speckle size below the bosonic healing length has recently al- 
lowed to clearly observe real-space signatures of Anderson localisation on the den- 
sity profile of a highly dilute gas expanding in a one-dimensional speckle potential 
[27j . The second method, to which our recent investigations have been devoted, 
is specific of optical lattice systems, and it is based on superposing a dominant 
(primary) standing wave with a secondary standing wave with an incommensurate 
period with respect to the primary one. In this way a quasi-periodic superlattice 
is created, whose physics interpolates between the one of a random system and the 
one of a perfectly periodic potential, as we will discuss in the next section. This 
approach has also allowed to recently observe clear signatures of Anderson localisa- 
tion in expanding one-dimensional atomic clouds with negligible interaction 58[ , as 
well as signatures of unconventional insulating behaviour in the strongly correlated 
regime |3a |. 



3.2 Ultracold bosonic atoms in disordered optical lattices 

Interacting bosons loaded in a deep optical lattice plus a random potential are 
described by a Bose-Hubbard Hamiltonian [45[ , of the form 

(io) h = -jJ2( b l b i + H - c -) + 1 E n ^ n - l ) + E( e > - 

(ij) » i 

where 6, is the boson annihilation operator, n, = b\bi is the boson number operator, 
J is the hopping integral, U is the on-site boson-boson repulsion, i runs over the 
sites of an optical lattice, and (ij) runs over the pairs of nearest neighbouring sites, 
e, is the on-site energy, accounting for the overall trapping potential (of the form 
Vt(i — *o) 2 ) and for any additional random or quasi-random potential. 

In absence of additional potentials, the Bose-Hubbard model at zero tempera- 
ture features a gapless superfluid (SF) phase with long-range phase coherence for 
large J/U; lowering this ratio leads the system to a quantum phase transition to- 
wards gapped Mott-insulating (MI) phases at integer fillings and with short-range 
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Figure 5 Schematic phase diagram of the Bose Hubbard model without disorder (left 
panel) and with disorder (right panel). MI=Mott insulator, SF=superfluid, BG=Bose 
glass, AG=Anderson glass. 




AG 



Figure 6 Phase diagrams for the Bose-Hubbard model in a Harper potential. Left 
panel: phase diagram for weak disorder A/U — 0.1 (in the inset the superfluid fraction 
and the compressibility gap are plotted as functions of t/U for fi/U = 0.25). Right panel: 
phase diagram for strong disorder, A/U = 0.5. 




coherence, (see Fig. [5]). Upon introducing a random on-site potential, namely tak- 
ing a as a random variable in the interval [— A/2, A/2], two new phases appear: 
for large J/U, an Anderson glass (AG) phase, in which weakly interacting particles 
are Anderson-localised by disorder, and their state connects continuously with that 
of the single-particle localisation problem; in the strongly interacting regime (low 
J/U ratio), a Bose glass (BG) phase, characterised by Anderson localisation of 
collective modes of the system, which can be essentially pictured as particle-hole 
excitations involving localised quasi-particle states. Both AG and BG are exotic 
insulators compared to the MI, in that they feature a gapless spectrum and a finite 
compressibility in absence of long-range coherence. The schematic phase diagrams 
of Fig. [5] can be obtained via mean-field theory [l(| and have been repeatedly con- 
firmed via numerically exact techniques, such as quantum Monte Carlo (QMC) [5rj| 
and density matrix renormalization group (DMRG) (5t1 |. 
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In the case of a quasi-periodic supcrlattice, the on-site potential takes the form 
of a Harper potential e« = A cos(2vrm + <j>) where < a < 1 is an irrational number 
and 4> an arbitrary phase factor. In the following we will specialise the discussion to 
one-dimensional systems, for which the problem of a single quantum particle in the 
Harper potential has been the subject of extensive investigations in the past (6^|. A 
one-dimensional particle in a purely random potential (as the one described above) 
is known to undergo Anderson localisation in all its energy eigenstates for any 
infinitesimal disorder strength, namely the Anderson transition occurs at A = 0. 
Moreover the single-particle density of states is continuous (on average over the 
disorder realizations). On the contrary a single particle in the Harper potential 
undergoes an Anderson transition from all extended to all localised eigenstates for 
a finite disorder, A = 2 J (26[. The density of states is not continuous (even upon 
averaging over the arbitrary phase </>), and it is structured in bands both on the 
extended and on the localised side of the transition. 

Our investigations have been focused on effect of strong interactions on the lo- 
calisation phenomena of bosons in a Harper potential, following a setup which has 
been recently realised in optical lattice experiments [36j. Making use of numerical 
exact diagonalisation, QMC and DMRG methods, as well as of analytical bosoniza- 
tion techniques, we have addressed the fundamental questions on the topology of 



the phase diagram in presence of a quasi-random potential ll|; 1 1 51 ] . Fig [6] shows 
the phase diagram of the Hubbard model in the Harper potential for two values of 
the disorder potential A/U = 0.1 (obtained via DMRG) and A/U = 0.5 (obtained 
via QMC), and for a = 830/1076, which is chosen as the ratio of the wavelengths 
used in the experimental realization, Rcf. [36| . In the case of weaker disorder 
( A/U = 0.1) one observes the appearance of a Bose-glass phase with vanishing 
supcrfluid density and finite compressibility at the interface between the conven- 
tional Mott insulating and superfluid phase. Yet, at variance with the case of a 
truly random potential 5(| 53] , a direct transition is possible between SF and MI 
close to the tip of the Mott lobe (as shown in the inset of Fig. [5]). This fact is 
intimately related with the existence of a finite localisation threshold associated 
with the quasiperiodic potential: analogously to what happens at the single par- 
ticle level, when the intersite hopping exceeds a critical value the quasiperiodic 
potential has no localising effect on the collective modes of the system, so that the 
Bose-glass phase is not observed, and a conventional phase transition between MI 
and SF can take place with the onset of gaplcss and propagating collective modes. 
The single-particle localisation threshold A = 2 J is attained for J/U = 0.2, which 
is already to the left of the Mott-lobe tip; on the other hand interactions play the 
role of screening the quasi-periodic potential, so that even a smaller J is expected 
to mark the disappearance of the Bose glass. This is quantitatively consistent with 
our data, where the BG phase seems to disappear for J/U « 0.15. 

On the opposite side of the spectrum, a very strong quasi-periodic potential 
has the effect of fully destabilising the MI phase of the homogeneous system at 
integer fillings; indeed for A/U = 0.5 the amplitude of the fluctuations of the 
Harper potential (~ U) overcome the gap characteristic of the MI phase, which 
is upper-bounded by U. Hence the MI phase disappears from the phase diagram, 
substituted by novel insulating phases stabilised by the quasi-periodic potential. 
Such phases reflect the intermediate nature of the quasi-periodic potential: in fact 
gapless and incoherent BG regions alternate with gapped incommensurate-band- 
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insulator (IBI) regions, which are the quasi-periodic analog of conventional band 
insulators appearing in interacting bosonic models in a commensurate superlattice 
(60| . In the limit of a deep quasiperiodic potential, the IBI phase can be associated 
to the opening of a finite gap for addition of a particle in any of the potential 
wells of the Harper potential: in this case an incommensurate " shell" of localised 
state is filled, and to start filling the next shell a gap needs to be overcome. This 
phenomenon does not occur in the case of a truly random potential, because the 
potential wells are not bounded in size, and hence the "charging energy" of each 
well can be arbitrarily small, making the system gapless at all fillings. 



3.3 Detecting disorder-induced phases 

The above results show that a very rich showcase of quantum phases can be 
implemented by experiments with strongly interacting bosons in incommensurate 
optical superlattices. Nonetheless the fundamental question remains on how one 
can detect experimentally the occurrence of such phases and all the details of com- 
plicated phase diagrams. This question is ultimately related with the potential use 
of these physical systems as analog quantum simulators [39j], capable of providing 
the answer to open theoretical questions concerning models which are implemented 
literally in the experiment. 

In the specific case of novel insulators induced by quasi-disordered potentials, 
our recent activity has been focused on how to detect their most striking feature, 



namely the simultaneous absence of a gap and of long-range phase coherence |12[ . In 
particular the gapless excitations of a BG are particle excitations, hole excitations, 
and particle-hole ones, and they are probed directly by the compressibility, which is 
indeed finite in the BG phase. The trapping potential confining the atoms offers the 
possibility of probing the compressibility by studying the response of the system to 
the increase of the trapping frequency (" trap squeezing" ) . Indeed such an operation 
tries to push particles from the trap wings towards the centre, where new particles 
are accepted only if the local state of the system has a vanishing gap to particle 
addition. Hence measuring a finite response of the average density in the centre of 
the trap upon increasing the trap frequency amounts to probing the realization of 
a locally compressible phase in the trap centre. 

Knowing that this phase actually corresponds to a BG requires the additional 
information on short-range coherence. The phase coherence of trapped cold atoms 
can be measured with time-of-flight measurements, which reveal the momentum 
distribution and in particular the zero- momentum component (coherent fraction): 

(11) n*=o = ^I>fo>' 

ij 

where N is the total particle number. In particular, if a particle is added to 
a BG phase, it will occupy a Anderson localised quasi-particle state, and hence it 
will not increase the coherent fraction of the system. Consequently a BG phase 
realised in the centre of the trap should be identified with a finite response of the 
central density to trap squeezing and on the simultaneous absence of response in the 
coherent fraction. This is indeed what it is observed in QMC simulations on trapped 
bosons subject to a variable trapping potential, as shown in Fig. [JJ The realization 
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Figure 7 Trap-squeezing detection of Bose glass phases. Upper panel: average central 
density and coherent fraction as a function of the chemical potential (controlled by the 
trapping potential [13]) for a system of iV = 100 trapped bosons in an incommensurate 
superlattice with A = U/2 = 10J and ot — 830/1076; boxed regions show BG behaviour, 
namely finite compressibility without any significant increase in the coherence. Lower 
panel: trap-squeezing-induced migration of a particle from a localised state on the trap 
wings from a localised state in the trap centre. The shown data refer to a variation of the 
trapping potential Vt from 0.013J to 0.014J. 
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of such a behaviour occurs for values of the chemical potential which correspond 
to BG behaviour of the homogeneous system, so that this type of measurement 
allows to reconstruct the phase diagram of the system in absence of the trap. The 
major requirement of this proposal is the measurement of the average density over 
a region of a linear size of 5 — 10/im, corresponding to 10-20 lattice sites. This is 
indeed possible with current large- aperture optics [631 ]. 
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